Bloch oscillations" in the Mott-insulator regime 
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We study the dynamical response of cold interacting atoms in the Mott insulator phase to a static 
force. As shown in the experiment by M. Greiner et. al., Nature 415, 39 (2002), this response 
has resonant character, with the main resonance defined by coincidence of Stark energy and on-site 
interaction energy. We analyse the dynamics of the atomic momentum distribution, which is the 
quantity measured in the experiment, for near resonant forcing. The momentum distribution is 
shown to develop a recurring interference pattern, with a recurrence time which we define in the 
paper. 

PACS numbers: PACS: 32.80.Pj, 03.65.-w, 03.75.Nt, 71.35.Lk 
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1. Introduction. Recently much attention has been 
paid to the properties of Bose-Einstein condensates of 
cold atoms loaded into optical lattices. In particular, the 
experimental observation of the superfluid (SF) to Mott 
insulator (MI) transition has caused a great deal of ex- 
citement in the field. Note, that besides demonstrating 
the SF-MI quantum phase transition, the same experi- 
ment also rose the problem of the system's response to 
a static force (used in the experiment to probe the sys- 
tem). Obviously, the response depends on whether the 
atoms are initially prepared in the SF state or the MI 
state. The former case was investigated theoretically in 
recent papers 

Him (see also related studies |lllH). 
It was found that, similar to the case of non-interacting 
atoms, the static force induces Bloch oscillations of the 
atoms which, however, may be affected rather dramat- 
ically by the presence of atom-atom interactions. The 
latter case of MI initial state was analysed in Ref. 0| , 
and is also the subject of the present Brief Report. In 
particular, we address the evolution of the atomic mo- 
mentum distribution not discussed so far. We show that, 
in formal analogy with usual Bloch oscillations, a static 
force causes oscillations of the atomic momentum, how- 
ever, with a different characteristic frequency. 

2. The model and numerical approach. Like in 
our earlier studies Q, |M Q > we model cold atoms loaded 
into an optical lattice by the Bose-Hubbard Hamiltonian, 
with an additional Stark term: 



cii -\- h.c. 
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In Eq. Q J is the hopping matrix element, W the on-site 
interaction energy, d the lattice period, and F the mag- 
nitude of the static force. Throughout the paper we con- 
sider a one-dimensional lattice and assume, for simplicity, 
that the filling factor (number of atoms per lattice site) 
equals unity. Then, the assumption of a Mott-insulator 




FIG. 1: Schematic presentation of the dipole state. Provided 
the Stark energy (mismatch between the 'bold levels' in the 
figure) is equal to the interaction energy (distance between 
the 'bold' and 'thick' levels) the atom may resonantly tunnel 
in a neighbouring well, thus creating particle-hole excitations 
of the MI state. 



initial state implies J < 0.1711^. In what follows, how- 
ever, we shall be mainly interested in the limiting case 
J <^W. Under this condition, the excitation of the sys- 
tem is only possible if the Stark energy is approximately 
equal to the interaction energy. Indeed, for dF = W 
the atoms may resonantly tunnel in the neighbouring 
well, thus forming 'dipole' (in the terminology of Ref. 0) 
states (see Fig. 0. 

The first step of our analysis is to identify the resonant 
subspace in the system's Hilbert space (spanned by Fock 
states), i.e. the manifold of states resonantly coupled to 
the MI state. For example, for a finite lattice with L = 8 
and F > 0, the MI state |11111111) is coupled to the 
one-dipole states |20111111), 1 12011111), etc., which are 
coupled to two-dipole states |20201111), |20120111), etc., 
which in turn are coupled to three-dipole states, and so 
on. If J ^ W, one actually can neglect the other (non- 
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FIG. 2: Upper panel: Dynamics of the normalised mean mo- 
mentum {p{t) ^ p{t)/JN) ior L = N = 9., W = 0.0324, 
dF — W, and J — 0.0019. The solid line shows the exact 
dynamics, dashed line corresponds to the resonant approxi- 
mation. Time axis is scaled with respect to Tj = 2nfi/ J, 
which is the characteristic time-scale of the system. Lower 
panel: Dynamics of the mean momentum calculated within 
the resonant approximation for different system sizes - L = 8 
(dot-dashed line), L = 12 (dashed line), and L = 16 (solid 
line). It is seen that, for the considered time interval t < Tj, 
convergence towards the thermodynamic limit is reached al- 
ready for L = 12. 

resonant) states when considering the excitation of the 
system. The validity of this resonant approximation is 
illustrated in the upper panel of Fig. 12 where we com- 
pare the time evolution of the mean atomic momentum, 
calculated by using the complete basis (for chosen L = 8 
the total dimension of Hilbert space isAf= 6435), and re- 
stricted to the resonant manifold respectively {Mr — 47) . 
It is seen that the resonant approximation works pretty 
well already for J/W = 0.05. It is also worth noting 
that, in the resonant approximation and after scaling 
time t — !■ Jt/h, the only relevant parameter of the system 
is the dimensionless detuning 

X={W-dF)/J . (2) 

Let us briefly comment on our choice of periodic 
(cyclic) boundary conditions used throughout this paper. 
These are imposed on after a gauge transformation 
which leads to the time-dependent Hamiltonian 

H(t) = -i (^e--* ^li^i + '^•^•) + Y E ^'("'-1) ' 

(3) 

with lob — dF/h the Bloch frequency. Note that for 
periodic boundary conditions the quasimomentum is a 
conserved quantity and, hence, the dipole states can be 
excited only in coherent superpositions, with the same 



quasimomentum k = as for the initial MI state. In 
particular, for onc-dipole excitations this constraint de- 
fines the state 

1 ^ ^ 

|i?(l)) = ^^5'|2011...11), (4) 

1=1 

where S denotes cyclic permutations: 5*12011... 11) = 
1 1201... 11). It is worth noting that for two-dipole 
(three-dipole, etc.) excitations there are many different 
states \D{2)), not related to each other by cyclic permu- 
tation. In what follows, we refer to the states \D{m)) as 
the translationally invariant dipole states. 

We conclude this section by a remark on the thermo- 
dynamic limit L oo. Obviously, the dynamics of a 
system of finite size differs from the one of an infinite 
system. However, this difference emerges only after a 
finite 'correspondence' time. This is illustrated in the 
lower panel of Fig. |2 which shows the mean momentum 
p{t) for different lattice size L = 8, 12, 16. It is seen that 
an increase of the system size above i = 12 does not 
change the result and, hence, for the considered time in- 
terval t < Tj convergence of thermodynamic limit has 
been reached. 

3. Results of numerical simulations. This sec- 
tion reports the results of numerical simulations the sys- 
tem dynamics obtained within resonant approximation. 
In our numerical simulations we followed the scheme 
of present-days laboratory experiments, where one mea- 
sures the momentum distribution of the atoms by using 
the 'free-flight' technique. Precisely, after preparation of 
the initial state (cooling stage), the atoms are subject to 
a static force for a given time interval t (evolution stage). 
Then the static field, as well as the optical potential, is 
abruptly switched off and the atoms move in free space. 
Finally, the spatial distribution of the atoms is measured, 
which carries information about the momentum distri- 
bution at the end of the evolution stage. Repeating the 
experiment for different time intervals t, one recovers the 
time-evolution of the momentum distribution P{p, t). 

Figure |3| shows the time-evolution of the atomic mo- 
mentum distribution for an optical potential depth v 
equal to 22 recoil energies, and a static force strength 
F corresponding to A = 2.084. Note that the v uniquely 
defines the hopping matrix element J (J = 0.0038 re- 
coil energies, for v = 22) and, hence, the tunnelling time 
Tj = 2'Kh/ J. The amplitude v also defines explicite form 
of the Wannier states ipi{x) = "00(2; — dl) and, thus, the 
initial distribution Pq{p) — P{p,t = 0) of the atomic 
momenta. Indeed, since the initial state is MI state, 
|5'(t = 0)) = l-D(O)) = |11 . . . 11), the momentum distri- 
bution at t = is simply the squared Fourier transform 
of the Wannier state, as can be easily derived from the 
definition of the one-particle density matrix p^: 

p{x,x') = Y^PLm{t)'tpl{x)'lpm{x') (5) 
l.m 

(t) = i^mjdml^it)) . (6) 
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FIG. 3: Static force induced dynamics of the atomic momen- 
tum distribution for cold atoms in a one-dimensional optical 
lattice. The depth of optical the potential is 22 recoil ener- 
gies, the dimensionless detuning A — 2.084. A periodic change 
of the distribution is clearly observed. Feeble oscillations of 
the distribution are an artifact, due to the finite system size 
(L = 16). 

As time evolves, the momentum distribution repeatedly 
develops a fringe-like interference pattern. More for- 
mally, 

P{p,t) = Poik)fip,t) , (7) 

where f{p,t) ~ fip + PL,t) is a periodic function of 
the momentum with the period p^ = 2TTh/d defined 
by the inverse lattice period. For currently considered 
case A = 2.084, the function |7J is also (almost) periodic 
in time with the period Tx w O.SSTj. Note, however, 
that /(p, t) is in general not periodic (or quasiperiodic) 
in time, although some characteristic time scale prevails. 
This statement is illustrated by the upper panel of Fig. ^ 
where the temporal behaviour of the mean momentum 
is displayed for different values of the detuning A. One 
clearly observes the decaying oscillations of the mean mo- 
mentum, where both the period of oscillations and the 
decay rate increase as the detuning is decreased. 

4.Quasienergy spectrum approach. From the 
point of view of Quantum Optics, momentum oscilla- 
tions of cold atoms is a process of subsequent excitations 
of the translationally invariant dipole states \D{0)) ^ 
\D{1)) ^ {\D{2))} ^ .... As an overall characteristic 
of this process we consider the mean number of dipole 
states 

2 

^W--zEEHc™,,WI', (8) 

m— j 

where the |cm,j(i)P are occupation probabilities for dif- 
ferent dipole states. (Note that D{t) < 1, due to the 
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FIG. 4: Mean atomic momentum (upper panel) compared 
to the average number of dipole states (lower panel), as the 
functions of time, for different values of the detuning A = 
(solid line), 1.042 (dashed line), and 2.084 (dot-dashed hne). 

chosen normalisation.) The dynamics of D{t) for three 
different values of the detuning A is shown in the lower 
panel of Fig. 0] A strong correlation between the num- 
ber of excited dipole states and oscillations of the mean 
momentum is clearly observed. 

The above results of our numerical simulations can 
be qualitatively understood by analysing the quasienergy 
spectrum of the system. An explicite form of the effective 
Hamiltonian, who's eigenvalues define the quasienergy 
spectrum, immediately follows from ||2J) by employing 
the resonant approximation and is given in Ref. j^. 
Namely, using the notion of dipole creation operator 
dj — {l/\/2)al_-^ai ( which creates a 'hole' at site I and 
a 'quasiparticle' at site / — 1) the effective Hamiltonian 
reads 

Heff = A J2 d\di - ^ E {^l + ' 

with the constraint that neither there can be more than 
one dipole at one site {d\di < 1), nor two dipoles at 
neighbouring sites (hard core repulsion). As noticed in 
Ref. (9j, the eigenvalue problem for the Hamiltonian ^ 
can be mapped to the energy spectrum problem for a ID 
chain of interacting spins, where a number of analytical 
results is known. In particular, there is an Ising quantum 
critical point at Ac = —1.850, with qualitatively differ- 
ent ground state of the spin system below and above this 
value of A. It should be noted, however, that in the con- 
text of our present problem (dynamical response of the 
system), the eigenvalue problem for the effective Hamil- 
tonian ||2J) defines the quasienergy spectrum, where the 
notion of the ground state has no physical meaning. 

To understand the characteristic structure of the 
quasienergy spectrum it is convenient to discuss it first 
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FIG. 5; Evolution of the spectrum of effective Hamiltonian 
@ under variation of the dimensionless detuning A. (Lattice 
size L = 16.) 

for finite L. The result of a direct numerical diagonalisa- 
tion of Heff for L = 16 is presented in Fig. El This figure 
shows the position of the quasienergy levels, as a function 
of the dimensionless detuning A. To^ avoid possible con- 
fusion with a similar figure in Ref. 9] , we note that here 
only the states of the same translational symmetry as 
the MI state (i.e., the states with zero value of the quasi- 
momentum) are shown. It is also worth mentioning that 
the spectrum has reflection symmetry and, hence, when 
discussing the dynamics (rather than thermodynamics), 
only the case A > needs to be considered. 

Let us discuss the quasienergy spectrum in more de- 
tail. It is convenient to start with large positive A. For 
a large A the spectrum consists of separate levels (or 
bunches of levels), which in the formal limit A — > oo 
can be associated with the dipole states (or family of 
dipole states) \D{m)) with given m {m < L/2). The 



lowest level in Fig. (Sjis obviously the MI state, the level 
above the one-dipole state \D{1)), followed by the fam- 
ily of states \D{2)), etc. (For negative A the situation 
is reversed and the MI state \D{0)) is associated with 
the most upper level.) The key feature of the spectrum 
is the finite gap A between the quasienergy level \D{0)) 
and the rest of the spectrum, exiting for arbitrary values 
of A. It is precisely this gap, what defines the characteris- 
tic frequency of atomic oscillations seen in Fig.0] In the 
thermodynamic limit we have A 1.43 for A = and 
A |A| for |A| oo ^Tj. Let us also note that in the 
thermodynamic limit and for A « the remainder of the 
quasienergy spectrum is continuous and gapless. This ex- 
plains an irreversible decay of oscillations, although the 
decay rate (and its A-dependence) remains an open prob- 
lem. 



5. Conclusion. We considered the response of the 
Mott-insulator phase of cold atoms in an optical lattice to 
a 'resonant' static force. Here the term 'resonant' means 
that the Stark energy dF approximately coincides with 
on-site interaction energy W. Under this condition, the 
atoms can tunnel in the neighbouring wells of the op- 
tical potential, thus creating particle-hole excitations of 
the Mott-insulator state (the 'dipole states'). This pro- 
cess directly affects the atomic momentum distribution, 
which is usually measured in laboratory experiments. 
Namely, the momentum distribution repeatedly devel- 
ops an interference pattern with a characteristic period 
Tx- This period is uniquely defined by the tunnelling 
time Tj = ^ixfij J (J is the hopping matrix element) and 
the energy gap between two lowest quasienergy levels, 
which, in turn, is a unique function of the dimensionless 
detuning A = (W^ — dF)/J. In particular, for A = 0, 
one has T\ w O.GTTj, what, for e.g., for sodium atoms 
(Eu/h — 271 X 8.9 kHz) in optical potential of 22 re- 
coil energies [J /Eji — 0.0038) corresponds to approx. 20 
milliseconds. 
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